On the coarse-grained density and compressibility of a non-ideal crystal 
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The isothermal compressibility of a general crystal is derived and discussed within (classical) 
density functional theory. Starting from the microscopic particle density, we carefully coarse grain 
to obtain the thermodynamic compressibility and the long wavelength limits of the correlation 
functions of elasticity theory. We explicitly show that the isothermal compressibility is not the 
wavevector to zero limit of the (total) density correlation function. The latter also cannot be 
obtained from the static structure factor measured in a scattering experiment. 
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I. INTRODUCTION 

In crystals, where translational invariance is sponta- 
neously broken, strain enters as additional thermody- 
namic variable in the free energy. It describes the distor- 
tion of the solid and thus the trace of the strain tensor 
is connected to the change in density. In ideal crystals, 
where all atoms can be unambiguously assigned to lattice 
sites and all lattice sites are occupied, both are actually 
equivalent. In real crystals, point defects like intersti- 
tials and vacancies are present and density fluctuations 
also change by defect diffusion^. The particle number 
in a unit cell of the crystal can change by deformation of 
the cell (captured in the strain) or by motion of point de- 
fects through the cell. The presence of defects opens the 
question how density and strain fluctuations are defined 
in real crystals. Here, no one-to-one mapping of atoms 
to lattice positions is possible. Therefore, the displace- 
ment vector field, whose symmetrized (in linear approx- 
imation) gradient gives the strain, cannot be obtained 
from the individual displacements of the atoms from their 
lattice positions. Only recently microscopic definitions of 
strain and density fluctuations in real crystals were de- 
rived from the statistical mechanical description of real 
crystals, overcoming this difficulty^. The latter work fol- 
lowed an earlier suggestion by Szamel and Ernst^. 

An intriguing finding of the microscopic approach con- 
cerns the coarse-grained density field 6n(r, t) which en- 
ters into the free energy or, consecutively, the theory 
of crystal elasticity^. Even for arbitrarily large wave- 
lengths, particle density fluctuations with wavevectors 
close to all (finite) reciprocal lattice vectors contribute 
to the coarse-grained density field. This arises because 
at the macroscopic level strain fluctuations and defect 
density independently cause changes in the total parti- 
cle density, and strain fluctuations require contributions 
from wavevectors around finite reciprocal lattice vectors 
in order to be determined in Fourier space. In this contri- 
bution, we discuss this at first surprising finding within 
the framework of density functional theory. 

Based on the microscopic definition of the coarse 
grained variables of elasticity theory, we can ad- 
dress another intriguing question, originally raised by 



Stillinger— ~— . Namely, whether the structure factor in 
the limit of wavevector to zero is analytic and approaches 
the compressibility? We find that due to the long-ranged 
displacement correlations, the small wavevector limit in 
the correlation function of the coarse-grained density field 
actually is non-analytic and depends on the direction rel- 
ative to the crystal lattice. We derive these results from 
density functional theory and can thus put them on a firm 
microscopic basis. Thus, we generalize results obtained 
within the harmonic crystal approximation^. From stud- 
ies on two-dimensional crystals it is known that defects 
enter the expression for the isothermal compressibility 
in a complicated fashion^. We generalize these results 
to crystals of arbitrary symmetry. Correcting the equa- 
tions (A8) in the appendix of Ref. Q, we also derive 
relations between fluctuation functions and thermody- 
namic derivatives. Based on these results the elastic 
constants of crystals with point defects could be mea- 
sured by microscopy techniques in colloidal crystals^ or 
by simulations^. 

The paper is organized as follows: in Sect. II we re- 
call definitions and results from Ref. [2[ and derive ex- 
pressions for the fluctuations of displacement and den- 
sity fields in terms of microscopic quantities. To facili- 
tate application of the expressions, in Sect. Ill we invert 
these relations considering two sets of independent fluc- 
tuations, coarse-grained density and displacement field 
or displacement field and defect density^. In Sect. IV 
and the Appendices we show that the small wavevector 
limit of the coarse-grained density fluctuation function 
does not agree with the isothermal compressibility k. In 
Sect. V the reason for the disagreement and (static) scat- 
tering experiments are discussed. We conclude that these 
experiments do not allow to measure the compressibility 
in a crystal, in contrast to the familiar situation in liquids 
and gases^. 
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II. COARSE-GRAINED FIELDS AND THEIR 
CORRELATION FUNCTIONS 

In crystals, according to the Goldstone theorem, the 
spontaneous breaking of translational symmetry induces 
long-ranged correlations. Specifically, the displacement 
vector field u(r, t) possesses correlations which decay 
like the inverse distance. While the familiar expres- 
sion u(r,i) = J2i u l (i)<5(r — R 4 ), with the displacement 
u' (t) = r ? (t) — R l of particle i from its lattice site R z , can 
be used in ideal crystals, in real crystals, in which defects 
are present, this expression is invalid. In order to find the 
microscopic definition of u(r, t) and of the other fields of 
elasticity theory, the following approach was followed. 



A. Microscopic starting point in density functional 
theory 

The microscopic approach to find the displacement 
field in a real crystal starts from the particle density 
field p(v, t) (considering, for simplicity, a one-component 
crystal of point particles interacting with a spherically 
symmetric pair-potential) 



A 1 



p(r,t) = £<S(r -!■*(*)) 



(1) 



where N is the number of particles in the volume V . 
Spatial Fourier transformation gives fluctuations close to 
vectors g of the reciprocal lattice 



<5p g (q,i) =p(g + q,t) 



where 



p(k,t) = J d\e 



n s VS t 



N 



q,0 , 



r p(r,t)=Xy 



•(() 



and 



1 1 N 



V 



(2) 



(3) 



Here, the general wave vector k was divided up into recip- 
rocal lattice vector g and wave vector q, which lies within 
the first Brillouin zone; ( ) brackets indicate canonical av- 
eraging at fixed temperature T. The n g are the Bragg- 
peak amplitudes (Dcbye- Waller factors) which serve as 
order parameters for the spontaneous breaking of homo- 
geneity (translational invariance). 

Classical density functional theory (DFT)i£~— can 
then be employed in order to derive the correlation func- 
tions of the microscopic density fluctuations in Eq. ^ . 
The fundamental Ornstcin-Zcrnikc relation provides a 
connection to the inverse density-density correlation ma- 
trix Jgg'(q) 

Vk B T5^> = 5> Pg (q) <W(q)> J g ' g "(q) . (5) 
g' 



The (infinite-dimensional) Hcrmitian matrix J gg ' is the 
double Fourier-transform of the direct correlation func- 
tion c(ri,r 2 ). 



(6) 



V 



-J d^ij rrt- 2 e ?;g ri e" 18 ' r2 e iq ri2 (^jrj ~ c(ri , r 2 ) 



The direct correlation function is one of the central quan- 
tities of DF T 15 i 17 and is obtained as second functional 
derivative of the excess free energy J- ex with respect to 
the average density profile, k B T c(ri,r2) — — feiElL 



8n(r 1 )Sn(r 2 ) ' 

Given an (approximate) expression for the free energy, 
J gg ' can thus be taken as known. It constitutes the only 
input for the ensuing theory. 

In Ref. @ , the dissipationless Zwanzig-Mori equations 
for (microscopic) density and momentum density fluctu- 
ations are given in terms of J gg / . They are solved by the 
following ansatz for the microscopic density fluctuations 



<5p s (q, t) 



-in s g a du a (q, t) - n g , (7) 



"o 



with Greek indices denoting spatial directions; repeated 
indices are summed over (Einstein summation conven- 
tion). Equation (0 expresses that for wave vectors q 
within the first Brillouin zone, the four macroscopic fields 
Sn(q, t) and <5u(q, t) determine the hydrodynamic con- 
tributions of the microscopic density field. Because equi- 
librium correlation functions of the microscopic density 
fluctuations on the left hand side of Eq. (0 can be cal- 
culated within DFT, this ansatz allows us to calculate 
correlation functions of the macroscopic fields. More- 
over, explicit expressions for the coarse grained density 
and displacement fields follow using the two summations 



n \ - 



(8a) 
(8b) 



and the relation |"g| 2 g = 0. For completeness, we 
list here the results for the macroscopic densities. The 
coarse grained density is 



(9) 



g 



where the average particle number density uq = N/V and 
a normalization constant A/o = X)gl n g| 2 appear. The 
displacement field is 

u a (q,t) = n * s .9/3 Spg(<l,t) , (10) 



with a normalization matrix Af a p = ^ g \ n g\ 2 9agp- Note 
that in equilibrium (u) = 0, allowing us to replace u with 
Su for notational uniformity. 
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Equations (|9]) and (| 10[) express the macroscopic fields 
in terms of microscopic particle density (JTJ). It is in- 
triguing that contributions from all finite lattice vectors 
g =/= are present in the coarse-grained density. Even in 
the limit of vanishing wave vector, q — > 0, it is not suf- 
ficient to measure particle density fluctuations close to 
the center of the first Brillouin zone, in order to deter- 
mine the thermodynamic density field in crystals. Fluc- 
tuations from the regions around all lattice vectors con- 
tribute and describe how macroscopically strain fluctua- 
tions and defect density independently cause changes in 
the total particle density. 



B. Correlation functions 



x (v Q n(n))(v^n(r 2 ) 

Kp-ysq-yqs + 0{q 4 ). 



(13c) 
(13f) 



The small wavevector limit and the index-symmetries 
= \i-iia and A Q ^ 7< 5 = X8 ai s = ^afsSy = are 
discussed in detail in Ref. [2[. Their connection to ther- 
modynamic derivatives will be recalled in Eq. (|27[) . The 
obtained set of equations is best given in proper matrix 
notation 



Vk B T5u=\ 



/ (Sn* Sn) 
n' 1 






(8u*Ju p ) J 



i k 



^*(q) 


^(q) N 


,Mq) 


A^ 7 (q) y 



kj 

(14) 



It is now conceptually straightforward albeit somewhat 
tedious to derive the correlation functions of the macro- 
scopic fields in terms of expressions containing the direct 
correlation function. Using Eq. (J7]), one gets 



(Sp g (q)Sp g ,{q)) = 

( i it *, \ r / \ \ ,5n*(q)5n(q) 

n s n sf [g a g (6u a {q)6up(q)) + ( 5 ) 

\ n 

*, ^n(q) .,5n*{q) 

n n 



(11) 



W a (Su a (q) ) - i( 5up(cD)g p , 

no no / 



Inserting this into Eq. §5§ and with the help of the two 
summations ([8]). one obtains the following set of equa- 
tions 



Vk B T = 


(Sn Sn) 

2 v (q) 
n o 


(—Sup)p /3 (q), 
n 


(12a) 


3 = 


(5n*5n) 

2 Mq) 


(—du s )X* 50 (q), 
n 


(12b) 


Q = 


{8u* a 5up)np(q) - 


(<5<-X(q), 
no 

(5ul-)^(q). 

no 1 


(12c) 


Vk B T5 a ^ = 




(12d) 



The g-dependent constants of elasticity v(q), /i Q (q), and 
A Q ^(q) are defined in terms of integrals containing the 
crystal direct correlation function. 



Kq) = ^y~ J ddfl /^2n(r 1 )n(r 2 )e- it »' ri2 
[ , ; - c(ri,r 2 ) 



(13a) 
(13b) 



n(rx) 
«z/ + £% 2 ), 

M<0 = ^V-J ddri J d<V 2c(r 1 ,r 2 )(l 

x n(n)V Q n(r 2 ) (13c) 
^iu a0 q p + O(q 2 ), (13d) 
k R T 



A^(q) 



V 



dtn J d ( V 2 c(r 1 ,r 2 )(l 



with Latin indices i = 0, a. Obviously, the matrix of cor- 
relation functions of the macroscopic variables are given 
by the inverse of the matrix of the constants of elasticity 



^ /Sn*Sn\ 






(5u* a 5up) 



-Vk B T 



' Kq) 




v M«(q) 


A a/ 3(q) 



The single matrix elements are^ 



(15) 



[ 5 ^)=vk BT ( 1 - + ^ 

ni V v v 



X a f 



v 



(16a) 

= Vk B T (y - nl Xllftpj 1 = Vk B TK~\ 

(16b) 

(Su* a Sup) = Vk B T (\ aj3 - p a v~ x p*^) 1 = Vk B TH~p 

(16c) 

-l ,,* \-i 



= Vk B T 

(—Sup) = Vk B T 
n 

= Vk B T 

nn 

-(5u* a — )=Vk B T 



no 



Vk B T 



Aq/3 + ^a^-t K 



(16d) 
(16e) 
(16f) 

(16g) 
(16h) 



The second line of each expression is a non-trivial alter- 
native, which is here given for completeness. 

We thus reached our first aim of expressing the corre- 
lation functions of the macroscopic variables, total den- 
sity and displacement vector field, in terms of integrals 
containing the functional derivative of the free energy 
with respect to density, i.e. the direct correlation func- 
tion. Let us note in passing that translational and ro- 
tational symmetry^ yields the expected ^-divergences or 
-dependences of the correlations in the limit q —> 0, es- 
pecially (Su^Siin) ex l/q 2 follows from A Q( g(q) oc q 2 . 
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III. INVERSE RELATIONS 

While equations (fT6|) predict the fluctuations of the 
macroscopic coarse-grained density and displacement 
field based on the constants of elasticity obtained from 
the direct correlation function and thus the free energy, 
experimentally, the inverse relations are of interest: ex- 
pressing the constants of elasticity of the crystal in terms 
of measurable correlation functions. Two different sets 
can be used in experiment. Either displacement field 
and coarse-grained density fluctuations can be recorded, 



or displacement field and defect density. For reference, 
we provide the inverse relations for both cases in this 
section. 



A. Including coarse-grained density 

The matrix equation (|f4j) can be inverted in order to 
find the elastic functions z/(q), /i Q (q), and A Q( g(q) in 
terms of measurable fluctuation functions. Straightfor- 
wardly one finds: 



Kq) 

Vk B T 



Aqff(q) 
Vk B T 



Ma(q) 

Vk B T 



.6n*6n, , .611*611. 1 ,6n* . 
(— — > _1 + <— — > <— 5u «) ( Su J u p) 



n 



„ Sn. ,Sn*Sn. _i ,Sn* . . 
(Su a —){ — — } ( Sup) 



n 



n 



Sn*Si 



> - (^-6u Q )(6u*Ju p )- 1 (6 U * s — )) 1 = R-\ 
lln H nn I 



n 



"0 



c-l 
D a/3 



{6u* a 6up) 1 + (6m*<5u 7 ) 1 (6u* — )R 1 ( 6us)(6u* s 6up) 1 



n 



no 



c-i/x * Sn>. / 6n*5n 1 

(8u* a 6u p )- 1 (6u;—)R- 1 . 

no 



1 „ Sn. ,6n*6n._, . 
{5u* p —){ — ) 1 (17a) 



r 



n 



(17b) 
(17c) 
(17d) 
(17e) 
(17f) 



These relations enable one to determine the functions 
A 0| g(q), /i Q (q), and f(q) from measurements of fluctua- 
tion functions. 



B. Including defect density 

Although the relation between the constants of elas- 
ticity and the fluctuations of the coarse-grained fields is 
complete, it is worthwhile to consider a second set of vari- 
ables. So far the displacement field u a appeared in two 
different ways. It contributes to the coarse-grained den- 
sity, but it also appears as broken symmetry variable. In 
this section we introduce the point defect density c in lieu 
of the coarse-grained density, and keep the displacement 
field. 

The correlation functions of the coarse-grained density 
and displacement field are easily transformed into corre- 
lations between the fluctuations of displacement field and 
the point defect density using the definition 



<5c(q, t) = —6n(q,t) — inoq a 6u a (q, t) 



(18) 



The set of variables {<5c(q), 6u a (q)} may be more appro- 
priate to describe an experiment when few defects are 
present and <5c(q, t) can be measured easily. Thus (jT4")) is 



transformed into 
VknTSu = 



(6c* 5c) 






(8u* a 6u a ) 



ik 



i/*(q) 


noV 5 (q) 


noV„*(q) 


A; 5 (q) 



(19) 



The combination of the constants of elasticity appearing 
here is directly connected to the hydrodynamic equation 
of the momentum density expressed in terms of point 
defect density and displacement field as hydrodynamic 
variables^. There the time derivative of the momentum 
density couples to the displacement field via (the negative 
of) 

A Ql s(q) = A a(3 (q) - iq a fip(q) + i^* a (q)qp + q a v(q)qp. 

(20) 

The coupling to the point defect density is given by (the 
negative of) 



V a (q) = — (V*(q) - iq a v{<\) 
n \ 



(21) 



The individual matrix elements of the correlation func- 
tions in terms of f(q), V^(q), and A 0( g(q), and the limit 
q — > 0, may be determined according to the steps in the 
previous paragraphs. In the same way, the transformed 



5 



constants of elasticity can be connected to measurable 
fluctuation functions. The comparison of the matrices 
in Eqs. (fT4|) and (fT9|) indicates the straightforward re- 
placements in Eqs. (|16[) and (JTTJ) . The dynamical matrix 
A a p(q) determines the wave equation of the momentum 
density, and its eigenvalues give the (acoustic) phoiion 
dispersion relations. 



IV. ISOTHERMAL COMPRESSIBILITY OF 
CRYSTALS 

In this section the compressibility is derived from a 
thermodynamic consideration, details are given in the 
Appendices. The situation described is one in which no 
pre-stress is applied to the crystal in equilibrium. 



A. Connection to elastic coefficients 

Starting with the definition of the compressibility of a 
fluid 



]_dV 
V~dp 



N 



(22) 



the question arises how this has to be generalized to de- 
scribe the additional degrees of freedom of a periodic 
crystal. The free energy of a crystal 



dF = -pdV + fidN + hapdU, 



(23) 



includes a term with a stress tensor h a /3 at constant 
volume V and particle number N times an extensive 
strain tensor U a p = Vu a p (The work done is SW = 
J h a p5u a pdV with the symmetrized linear strain tensor 
u a /3 = \(S a u [i + V/3 u a))- The compressibility for a pe- 
riodic crystal is understood as the derivative at constant 
stress tensor 



]_dV 
V~dp 



N,h af! 



Tig dfi 



1 dn 




1 dn 




( dh-fS 




1 dhys 


Hg dfi 




nl du a/ 3 




\du a p 


J 


dfi 



K0 

(24) 

using Maxwell and Gibbs-Duhem relations described in 
Appendix [B] There this result is further manipulated to 



n 



dh 7 $ 



dh 



dn 



7<5 



dn 



/dfx 



(25) 



which is one of our main results. This expression for the 
isothermal compressibility of a general crystal generalizes 
results obtained for high symmetry crystals^, and, to 
our surprise, appears novel. Now, the expressions of the 



clastic constants as thermodynamic derivatives^ can be 
entered. They follow from the Gibbs fundamental form 
of the free energy density / = F/V given by 



df = udn + h a pdu a j3 , 



(26) 



and take the microscopic expressions already introduced 
in Eq. (CG 



9V 

dn 2 u aft 

d 2 f 
dndu a p 

d 2 f 



dfi 

dn 



diiapdu-fS n du 7 s 



du a p 
dh a p 



v/nl, 
dh aj3 



(27a) 

-Hap/no, (27b) 



dn 

(27c) 



See Appendix[A]for an alternative derivation of ([27| from 
the second order variation of the free energy in DFT as- 
suming smoothly varying density profiles. It supplements 
the derivation in Ref. [2[ (recalled in Eq. (jTB")) ). where 
the equivalence of the hydrodynamic equations was used. 
The compressibility is thus 



v Va/3 



C 



fl^gl 



(28) 



Neglecting the coupling tensor u a p = the second term 
vanishes and the compressibility k is given by kT 1 = 

, which plays the role of the inverse bulk mod- 



l dn 



ulus at constant strain. 



B. Including defect density 

If one considers the set of independent variables with 
the defect density c instead of the coarse-grained density 
n with Eq. (|18[) simplifying to dn = —nodu aa — dc the 
manipulations leading to Eq. (|2"5"|) have to be adapted. 
The compressibility is given now in terms of derivatives 
of the stress tensor at constant defect density^ a a p and 
the chemical potential fx, which are functions of the strain 
tensor and the defect density (df = —fide + a a pdu a p) 







-i dn 




X 


u a p' Van 




du a /3 


71 





2 <9/J 

l °dc 
( da-ys 



dfi 



du 



a/3 



du a 



da-ys 




(da 




1 dfi 




c dc 




[dc 




du a f3 


• 



X <5 7 5 + 



1 da^s 



no dc 



d/i. 
dc 



(29) 



The relevant thermodynamic derivatives are now given 
by (see Eq. (|A15[) for the free energy density) 



10 'd# 



,dfi 



(30a) 



6 



d 2 f du da a R 

n Q Q — = -no— — = n Q = vb a p + /i a p, 



dcduap 
d 2 f 



du af: 



dc 



(30b) 



du a pdu~fS 



(30c) 



An interesting limit is now the vanishing of the cou- 
pling between the defect density and the strain field, 
a a 1 — = 0. This yields two independent contributions 
to the compressibility 



(31) 



The first contribution v is due to the fluctuations of the 
defect density, whereas the second one (C^^g) -1 S a pSjs 
is due to independent fluctuations of the strain tensor, 
which agrees with the known result for a perfect crystal 
without external strain^. 



V. CONNECTION - SMALL WAVEVECTOR 
LIMIT 

So far we considered correlation functions and the 
isothermal compressibility of crystals. In this section we 
bridge the gap between the density correlation function 
and the compressibility, and point out the subtle differ- 
ence between the two expressions. In the second part of 
this chapter the so called generalized structure factor is 
discussed. 

In order to understand the connection to the compress- 
ibility, the g-dependence in the limit q — > of the corre- 
lation function of the coarse-grained density (|16a[) needs 
to be discussed in detail 



Sn*Sn 
Vk B Tnl 

A 

1 Mqflgft 
V V 



}=^ 1 (q) + ^ 1 (qK(q) 



x (Aa 7 (q) 1 (q)A*7(q) 



Mq> _1 (<i) (32) 



(33) 



Here we used the known small- wave vector expansions of 
the elastic coefficients, which were defined in Eqs. (|T3| . 
They follow from DFT relations expressing translational 
and rotational symmetry^. Noting that only the sym- 
metrized combinations in a O 7 and e (f> of the term 
in square brackets are relevant, and with the help of 
Eqs. ([27]) this expression can be further simplified to 



, 5n*5r 
'Vk B Tn z 



9^0 1 
2/ — ,," 



l°ac70 ~ lit 



1 n-fsqs 

V 

(34) 



This expression would agree with the thermodynamic 
one (p?5|) . if the factors of qpqs canceled q t q^,- That the 



limit q — > is not that simple can be seen even for highly 
symmetric crystals. For a cubic crystal, the correlation 
function yields different results in the small q limit (|34|) 
depending on the direction of q relative to the unit cell. 
And for the hypothetical model of an isotropic crystal, 
the small q limit (|34|) is direction independent, but dif- 
fers from the thermodynamic value from (|28p . The lat- 
ter simplified case, allows to identify the origin of the 
discrepancy and will be studied in detail in the next sec- 
tion. Appendix [D] considers the non-analytic nature of 
the small-q limit from the point of view of the harmonic 
free energy derived in Appendix [A] 



A. Perfect crystal embedded in a matrix 

To study the difference in more detail, it first 
simplification, more convenient to look at the simpler 
problem of a perfect crystal. In this section we also use 
the more familiar expressions of elasticity theory. The 
connection to the terms used so far is given in Appendix 
ICl For a perfect crystal the correlations of the displace- 
ment field is given by the (inverse) of the dynamical ma- 
trix D a p(q) 



(Su*Ju p ) = VksTD-Uq) 



(35) 



The coarse grained density for a perfect crystal is 
<5n(q, t) = —inoq a 5u a (q,t) and the dynamical matrix 
is related with the elastic constants^ via Z? Q7 (q) = 
Cap-ysqpqs- Thus for the coarse grained density corre- 
lation function we get 



, 6n*6n 
"Vk B Tnl 



) = QaDAicfiqp = q a (C a e^qeq^>) 1 qp (36) 



which shows the same problem in the limit q — > as 
arises in Eq. (|34p . when compared to the thermody- 
namic compressibility-^ n = (C~^ 7(5 )5 a/ 3(5 7( 5 (contraction 
of the inverse of the matrix of elastic constants). For an 
isotropic crystal the elastic tensor simplifies to the two 
Lame coefficients C Q( g 7( 5 = \5 a p5 1 s + u(5 a7 6p5 + SasS^). 
Thus, the compressibility is k^ 1 = A + whereas the 
correlation functions yields X + 2/i (which corresponds to 
the longitudinal speed of sound). 

To show the origin of this difference we consider an 
isotropic (ideal) solid for which the so called fundamen- 
tal solution of elasticity is known. Other symmetries with 
known solutions are hexagonal^ and pentagonal^. The 
cor resp onding problem in two dimensions can be found 
in |22|. We consider a three dimensional sphere with 
volume Vb embedded in a spherical matrix V of the 
same isotropic material. The radius Rb of the embedded 
sphere is increased Rb — > Rb + Ar and the surrounding 
matrix is compressed. To determine the displacement 
field and the elastic energy of such a deformation one 
has to solve the equation of elastostatic theory, which is 
the vanishing of the divergence of the stress tensor, or in 
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terms of displacement field 

V / 3C Qi a 7 aV 7 U5 = 0. 



(37) 



The solution is a sphere with increased volume Vb + AVg. 
The only non-vanishing displacement field is (homoge- 
neous dilatation) 



A/ 



Rb 



r < R B 
r > R B 



This yields for the total elastic energy 



E 



Vb/AVb^- 
2 V V B 



Vb 
V 



(38) 



(39) 



The first part is due to the stretched sphere and the sec- 
ond contribution is from the surrounding matrix. Thus, 
depending on the ratio the relevant combination of 
elastic constants changes from A + |/z for — > 1) to 
A + 2/j, for — > 0). In the limit of vanishing shear 
modulus n the difference vanishes. Thus, for a fluid it 
doesn't matter if one determines the volume fluctuations 
of a small sphere in surrounding fluid (of the same kind) 
or if one looks at the global fluctuations of the whole 
system. 

It is worthwhile to note that the same ratio between 
these two combinations of Lame coefficients appears in 
a related context. In Eshelby's study22 of an inclusion 
in a matrix of elastic material, the so called constrained 
strain it^ is given by the stress-free strain u^p 



\ _i_ 2 
C _ A + 3£_ T 



(40) 



An analog problem is a polar fluid in a dielectric 
mediu m 16 i 24 ' 25 . There the susceptibilities show a direc- 
tional dependence due to the dipolar interaction, and a 
different combination of dielectric constants is relevant 
depending on the surrounding medium. 



B. Generalized Structure Factor 

There is a further aspect which differs in the relation 
between the compressibility and the correlation of the 
density fluctuations of a fluid and a crystal. There is a 
difference if one looks at the elements of the generalized 
structure factor which contribute to the compressibility, 
i.e. whether those are different from S f g= o(q — > 0). 

We recall that the generalized structure factor is de- 
fined by 

S g (k) = I J <fV 2 (5p( ri )^(r 2 )) e- l s R e -' kAr 

(41) 



-<<yp( g /2 + k)*p(g/2-k)> , 



(42) 



(with R = (ri +r 2 )/2 and Ar = ri — r 2 ) and its So(g + 
q) element is measured in a scattering experiments^. 

In a liquid, where translational invariancc dictates 
that only So(k) is non- vanishing and isotropic, its 
connectio n 15 ^ 6 to the compressibility is well known 
So(q —> 0) — > rvoksTK. To convince oneself that such 
a connection does not hold in a crystal, the definition of 
the coarse-grained density Eq. ^ can be used to derive 



5n*5n 



(2tt 



72 "S' vV ( g ' + <l) 5 Pte + l)) 7 



4. 



A/? 



5ji g £)Sg(--g- q)£(g - g' - g)n* 

g:g' g 

N 2 Z^^'^g 



J g-g' 



g,g 



g + g 



q n. 



(43) 



where the left hand side becomes k for q to 
zero in the fluid case. Clearly, every element of 
(<5/9* (g' + q)5p(g + q)) is involved, not just the one with 
vanishing reciprocal lattice vector g = g' = 0. Also 
the correlation of coarse-grained density fluctuations is 
given by a combination of generalized structure factors 
S'g_g/(— (g + g')/2— q) in the limit q — > and not just by 
<S'g=o (q — > 0) as for a fluid. Although the possibility that 
the RHS of the last equation is indeed the compressibility 
cannot be ruled out, it seems rather unlikely. 



VI. CONCLUSIONS AND OUTLOOK 

We derived the thermodynamic expression for the 
isothermal compressibility k in a general crystal, and dis- 
cussed its connection to the limit of vanishing wavevector 
in the density correlation function. The correlation func- 
tions were calculated within the framework of density 
functional theory, with appropriate coarse-graining, in 
order to allow for a finite density of (point) defects. Ex- 
plicit expressions for the coefficients in the phenomeno- 
logical free energy in terms of the direct correlation func- 
tion of density functional theory were obtained. The cor- 
relation function of the coarse-grained density field does 
not, in general, determine the compressibility. For the 
case of an ideal isotropic solid, we could identify the ori- 
gin of the discrepancy from a calculation in macroscopic 
elasticity theory. It arises from the long-ranged strain 
fluctuations which decay like 1/r 3 and thereby cause 
boundary effects to enter the elastic energy. While in sys- 
tems with spontaneously broken symmetry, anomalous 
longitudinal correlations exist in general^ (besides the 
familiar symmetry restoring fluctuations^), the present 
observation appears more related to long-ranged dipolar 
correlations in polar fluids^. There the dielectric tensor 
in response to the vacuum electric field depends on the 
shape of the material and on the boundary conditions. 
It can be connected to a well defined isotropic dielec- 
tric constant only via shape/ boundary-effect dependent 
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distribution functions. For crystals, it remains to un- 
derstand this relation for arbitrary symmetries, and to 
derive it from a discussion of the microscopic correlation 
functions. 
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Appendix A: Free energy 



n(r 1 )n(r2)fa(r 1 )fa(r 2 ) , . , ^ n(r 2 )Sn(r 2 ) 

2 Hay^l J* a^v**!/ 



(2.) 



n(ri)(5n(ri) 



»o 



u a (r 2 )V a n(v 2 ) , 



(4.) 



(3.) 



(A3) 



which is, in the following, term by term manipulated with 
the help of the LMB^W— equation 



The second order change in free energy AJ 7 due to 
a variation in the density distribution Sp(r) around the 
periodic crystalline equilibrium density i o 16 ' 17 ' 27 

P*F=\ J J ^1^(^-0(1!,^))^^)^), 

(Al) 

where c(ri , r 2 ) is the direct correlation function of a peri- 
odic crystal. Assuming a smooth density variation, which 
is the real space expression of (0, i.e. (with the nota- 
tional simplification i5u(r) = u(r) ) 



5p(r) = — u(r) • Vrt(r) + n(r) 



Sn(r) 
n 



(A2) 



one can evaluate the thermodynamic derivatives appear- 
ing in Eqs. ([27]). Thus 

5p{r 1 )5p{v 2 ) =UQ(ri)u / 3(r 2 )V a n(ri)V / 3n(r 2 ) 



Vqn(r) 
ra(r) 



= / dVc(r,r')V a n(r'). 



(A4) 



The first part is discussed explicitly for a crystal in 
Masters^ and is eq uivalent to the discussion of the sur- 
face tension in [31| . The other parts are transformed 
accordingly. 

One expands u^(r 2 ) around i"i, which is valid for a 
short range (in ri 2 ) direct correlation function. The 



zero order term of AJ 7 ^ \ vanishes, because of (|A4[) and 
the first order term does not contribute as c(ri,r 2 ) = 
c(r 2 , ri). Since the hydrodynamic variable u(r) is slowly 
varying, one obtains an expression which is quadratic in 
Vu(r) 



(i.) 



A ^(i.)= ^T" // ^1^2(^4 -c(r 1 ,r 2 )) Ua (r 1 )^(r 2 )V a n(r 1 )V^n(r 2 ) (A5) 



2 JJ -Vn(n) 
k B T 



' j J ^i^V 'gn(r 1 )c(r 1 ,r 2 )\/ 'ffn(r 2 )M a (r 1 )^ff(r 1 )-u^(r 1 ) + y 7 Uff(r 1 )r 12 ^^, -^V 7 V^ )3 (r 1 )r 12i7 r 12 ^^ 

—0 symmetry 

= - J d ( VA Q( 3 7 «5V 7 u Q (r)V5U / 3(r), (A6) 
J J d t Vid < V 2 V Q n(r 1 )c(r 1 ,r 2 )V ( gn(r 2 )r 12 , 7 ri 2 ^ (A7) 



, _ k B T 



The constant of elasticity X a p^s is the same one as was 
obtained in the framework of hydrodynamic equations of 
motion 2 -. For the second term of the free energy, Sn(r 2 ) 
is expanded around ri and, as hydrodynamic variable, 
assumed to be slowly varying 

AT {2 , = k -^J j rfS^ /" (ri )f (r2) 



x [n(ri)6(ri 2 )-n(ri)c(ri,r 2 )n(r 2 )], 

(A8) 
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With 



V 



J jdfirxdfrz («(ri)5(ri 2 )- n(ri)c(ri,r 2 )n(r 2 ) 

(A10) 

The third and fourth part yield with the same arguments 

5n(r) 



AT, 



(3.+4.) 



= - J Sr ^ 



-V f3 u a (r), 



(All) 



n a p : 



k B T 
V 



J d'Vi J d < V2«(ri)V Q n(r 2 )ri2^c(ri,r 2 ) 

(A12) 



So all of the constants of elasticity v, /i a p , and X a p-yS 
appear in the expression for the free energy 



AT 



\ [ ^K^) 2 + C ^ V/3Uc,(r)V5 ^ (r) 

(A13) 



S tl ( r 1 

V/3U a (r) 

n 

with the Voigt-symmetric elastic constants C™p s — 
A a7 /3<5 + Xp-yaS — Xap-ys- The replacement in the last 
step reflects the rotational symmetry as only symmet- 
ric combinations of strain enter into the elastic energy 
and is based upon the rotational analog of the LMBW 
equation^! 

rixV (1) lnn(ri) = J d'^c^x, r 2 ) (r 2 x V (2) n(r 2 ) 

(A14) 

Let us add that the above expression for the free en- 
ergy also determines the constant Cq = in Eqs. (89), 
(90), and (92) of @, which could not be determined in 
where the elastic constants were identified via the 
hydrodynamic equations. 

The free energy in terms of the defect density <5c(r) 
and the displacement field u(r) is 



AT 



L f / oc( V } \ ^ / \ 0C\ Y I 

- cfirv — — +2(i/6 a p+fj, a p)—^-Vpu a {T) 

A.I \ Tin / V / HQ 



6c{r)\ 2 

2 J ur "\T^ 

(A15) 

Appendix B: Isothermal compressibility 

As a consequence of (|23p a Gibbs-Duhem relation can 
be derived 

- Vdp + Ndfi + U aP dh a0 = 0, (Bl) 
which yields for isothermal change with dh a p = 

Ndfi = Vdp. (B2) 

With fixed volume 



A d ^ 

d(1= dN 



dp 

dN, and dp = — 

V,h al3 dN 



V.h c 



dN. (B3) 



This yields 
du 



N 



dN 



v,h aB dN 



dp 

v,h afl dn 
N dV 

V 2 dp N,h 



V,h a p 

-1 



= (no/s) . (B4) 



This is further manipulated with the Jacobian trick^ to 
yield Eq. 



1 dn 
Hq du 



dn 
1 dp 


dh-ys 


dn 


dh^s 
du 




"a/3 P 


du a $ 


n 2 dh l5 
du a p 





1 dn 




1 dn 




( dh 7 s 


.) 


1 9/l 7 5 


n\ du 




nl du a p 


A* 


\du a p 




dfi 



(B5) 

The expression in the denominator is understood as the 
inverse of a fourth rank tensor A Q/ 3 7 5, which is given by 



(B6) 



The first term of the compressibility is basically the only 
non- vanishing term in the fluid limit 



J_dn 

nl dfi 



dfi 

dn 



(B7) 



For the second term the chemical potential /i is expressed 
as a function of density and strain tensor u{n, u a p) 



du 

which yields 



dfi 

dn 



dn + 



dp 



dUafj 



du 



a/3 



(B8) 



dn 


- _(!*E 




-i du 




du a p 


u \dn 


u a/3 J 


du aj3 


n 



The last two terms are 



dh^s 




dh 7 s 


dh^s 


dn 


du af3 




du a p 


n dn 


u a p du a p 



dhjg 




dhjs 




(du 


y 


du 


"off 


dn 


"of) 


ydn 





-n u fiys- 



= nov V«/3- (B9) 



(BIO) 
(Bll) 
(B12) 
(B13) 



Thus the compressibility of a periodic crystal in terms of 
its elastic constants v, fj, a /s , and C™^ 7(5 is indeed Eq. 



K — v 1 + v 1 u a p 



U^V 



(B14) 
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Appendix C: Elasticity 



With the expressions of Sect. lIIIBl Eq. (|36|) reads 

(5u*Ju p ) = Vk B TA-^a). (CI) 

This follows from (fT9f with V a (q) = 0, or v8 a p + ii a (3 = 
in the low q- limit, for a perfect crystal without point 
defects and vanishing coupling. The correlations of the 
coarse-grained density (|35[) become 



, Sn*Sn 
'Vk B Tn2 



) = ?aA a g(q)g l s = g a (C Q 



_1 W. (C2) 



where we used the small q expansion of the constants of 
elasticity and took care of the proper symmetric combi- 
nation as discussed in @- The elastic constants C a p 7 5 
of elasticity theory correspond to C^^ s in Eq. (|30cl) . 
which, with the help of v8 a p + fi a p = for a perfect 
crystal without coupling, could be further simplified (to 
C a pf5 = — v5 a j35 1 s). 

The elastostatic theory is contained in the static limit 
of the hydrodynamic equations of motion, see Eqs. (87) 
in 0. Without point defects the only non- vanishing 
equation is (87c), which then reads 



(C3) 



But this is just the Fourier-transformed equation of elas- 
tostatics (13711. 



Appendix D: Equipartition principle 



The harmonic free energy Eq. (|A13|) from Appendix 
can be written in a more compact form with the help of a 
4 x 4-matrix of elastic coefficients introduced in Eq. (fT4"]) . 
Fourier-transformation leads to 



AT 



<?»*(q) 

, n 



d d q 
(2n) d 

, £<(q) 



(Dl) 



<Mq) 

no 

<5w 7 (q) 



The relation ([14]) , which follows from the ansatz ([7]) , thus 
is identified as statement of the equipartition theorem re- 
sulting from this quadratic free energy. It leads to the 
wavevector-dependent correlation functions of the coarse- 
grained fields for small q. 

Interestingly, also the thermodynamic fluctuation func- 
tions, which correspond to q = 0, can be obtained from 
the harmonic free energy Eq. (|A13[) . even though they 
were not the values reached by the limit q — > in the 
q-dependent correlation functions. Equation (|A13[) can 
be rewritten using the Voigt notation, which shall be 
given by latin indices 1 < i < 6, where Ui = u a p for 
(a, P) = {(1, 1); (2, 2); (3, 3); (2, 3); (1, 3); (1, 2)}. For spa- 
tially constant fluctuations indicated by subscript q = 0, 
one obtains using the thermodynamic derivatives from 
Eq. (USD 



A J" 



V 



Sn 
n ' 



(D2) 



Now the 7x7 matrix of covariances is given by 



(SnSn) 




= Vk R T 



v -fj,j 
-Mi C™ 



(D3) 



The correlation of the density fluctuations is thus 



,SnSn, 



= vk B r(- 
= vk R r(- 



v 

Ha/3 
V 



V I 



v r 

(D4) 



where in the last line the usual notation is used instead 
of the Voigt one. This equipartition result agrees with 
the compressibility in Eq. (|28|) . 

Reassuringly, it suffices to discuss the harmonic approx- 
imation to the free energy in order to observe the non- 
analytic structure of the correlation functions in a crystal. 
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